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Abstract 



Direct calculation of the one-loop contributions to the energy density of bosonic 
and super symmetric (j) 4 kinks exhibits: (1) Local mode regularization. Requiring the 
mode density in the kink and the trivial sectors to be equal at each point in space yields 
\ the anomalous part of the energy density. (2) Phase space factorization. A striking 

position-momentum factorization for reflectionless potentials gives the non-anomalous 
energy density a simple relation to that for the bound state. For the supersymmetric 
kink, our expression for the energy density (both the anomalous and non-anomalous 
parts) agrees with the published central charge density, whose anomalous part we also 
compute directly by point-splitting regularization. Finally we show that, for a scalar 
field with arbitrary scalar background potential in one space dimension, point-splitting 
regularization implies local mode regularization of the Casimir energy density. 
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Introduction. — Quantum corrections to solitons were of great interest in the 1970's and 
1980's, and again in the last few years, due to the present activity in quantum field theories 
with dualities between extended objects and pointlike objects. Dashen, Hasslacher, and 
Neveu pQ, in a 1974 article that has become a classic, computed the one-loop corrections to 
the mass of the bosonic kink in </> 4 field theory and to the bosonic soliton in sine-Gordon 
theory. For the latter, there exist exact analytical methods associated with the complete 
integrability of the system, authenticating the perturbative calculation. Our work here uses 
general principles but focuses on the kink, for which exact results are not available. Dashen 
et al. put the object (classical background field corresponding to kink or to s-G soliton) in a 
box of length L to discretize the continuous spectrum, and used mode number regularization 
(equal numbers of modes in the topological and trivial sectors, including the zero mode 
in the topological sector in this counting) for the ultraviolet divergences. They imposed 
periodic boundary conditions (PBC) on the meson field which describes the fluctuations 
around the trivial or topological vacuum solutions, and added a logarithmically divergent 
mass counterterm whose finite part was fixed by requiring absence of tadpoles in the trivial 
background. They found for the kink 

M (1) = E - E ~^i 0) + AM = -hm (A _ ^ < o , (1) 

where m is the renormalized mass of the meson in the trivial background. This result remains 
unchallenged. 

The supersymmetric (susy) CclSC, clS well as the general case including fermions, proved 
more difficult. The action and the kink solution read 



l/o ,\2 1-To , 1 dU- 



0kink(a;) = -^tanh^ , (2) 

where — \U 2 = — \(<t> 2 — /i 2 /A) 2 , the meson mass is m = //\/2, and c = 1 for supersymmetry. 
Dashen et al. did not explicitly compute the fermionic corrections to the soliton mass, stating 
"The actual computation of [the contribution to] [due to fermions] can be carried out 
along the lines of the Appendix. As the result is rather complicated and not particularly 
illuminating we will not give it here" (page 4137 of [Q). 

Several authors later computed M*- 1 -* for the susy kink and found different answers. It 
became clear that the method of Dashen et al. yielded results that depended on the boundary 
conditions for the fluctuating fields. In fact, repeating exactly the same steps for the susy 
kink as taken by Dashen et al. for the bosonic kink (using PBC also for the fermions, 
i[)±(—L/2) = ip±(L/2)), taking equal numbers of modes in all four sectors, including one 
term with uj ~ (due to a periodic solution with u 2 > 0) in the bosonic kink sector and 
one term with u = in the fermionic kink sector (explicitly, there are two real independent 
solutions with u — 0, one localized at the kink and one at the boundary, and the coefficient 
of each satisfies c 2 = | |2J gives = hm(j — j-) [3]. We now know that this is the 
correct answer to an inappropriate question, because it includes boundary energy. Schonfeld 
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[3] finessed the problem of a single kink with its sensitivity to boundary conditions by 
considering the kink-antikink system with PBC. Taking into account two terms with uo ~ 
in the bosonic kink-antikink sector (due to one periodic solution with u 2 < and another 
antiperiodic solution with u 2 > 0) and one term with uj = in the fermionic kink-antikink 
sector (due to two periodic solutions with uj = [2J, he obtained what we now know to 
be the correct answer. The problem of boundary contributions was circumvented by other 
methods in [3] and The fermionic contribution to is given by 

>0 , (3) 




and the total one-loop correction is thus 



MW =-?™ . (4) 



2tt 



With attention restricted to the kink alone, it was shown in [2] that one could eliminate 
boundary contributions from the fermionic part of the energy by averaging over quartets of 
boundary conditions for the fermionic fluctuations - periodic, antiperiodic, twisted periodic 
and twisted antiperiodic, where twisting means interchange of the upper and lower compo- 
nents of the fermion wave function ip± — > ip T 5j. This averaging is necessary to preserve 
certain discrete symmetries for fermions. The results in [2] give a complete, though intricate, 
way to calculate in terms of mode frequencies u n . 

In light of the complexities which boundary conditions generate for the problem including 
fermions, the most important advance since [T| was the approach of Shifman, Vainshtein, and 
Voloshin j.7 , who used higher space-derivative regularization (with factors (1 — d 2 /M 2 ) for 
the kinetic terms but not the interactions) to compute the central charge densities of the susy 
sine-Gordon soliton and kink. Their scheme is manifestly susy, canonical (no higher time 
derivatives), and independent of boundary conditions (because it yields a local density). 4 
They argued that the energy density is equal to the central charge density (because the 
difference is a susy commutator) and they computed the latter - including an anomaly 
recognizable as an M 2 /M 2 effect. They verified that the one-loop correction to the 
integrated central charge of the kink comes only from the anomaly and is equal to (j3J). The 
presence of a topological anomaly was first conjectured in j^J. 

One may compute the energy density for the bosonic sine-Gordon soliton by mapping 
the system onto another one which exhibits supersymmetry, and computing the central 
charge density for that fictitious system [7j. Quite possibly similar techniques would work 
for the bosonic kink. Our approach here is instead to attack the Casimir energy density 
directly, freeing the calculation from dependence on supersymmetry. In doing so, we have 
found it helpful to introduce a simple rule for regularization of energy densities which we 
propose as a fundamental principle. At this point, the primary evidence for the validity 
of the principle is the agreement between the energy density we compute by its use with 
the central charge density of |7,. Our proposed principle appears new in the literature, 

4 As [J] point out, Yamagishi |H] was perhaps the first to study local densities in this context. 
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and yet has roots in early quantum physics. Local mode regularization (lmr), or mode 
density regularization, is the local counterpart of the familiar (global) mode regularization 
or mode number regularization. The principle is that, when fluctuations are expanded into 
normalized modes 4> n (x) (i.e., / dxcf)* n (x)4> n (x) = 1 for all n), the cut-off local mode density 
Pn(x) = J2n=i^n( x )^n(x) should be background-independent, meaning the same in the 
trivial (p^) and kink (p) sectors. For the bosonic kink this implies that one must truncate 
the sums at different upper bounds corresponding to different frequency cutoffs 



Here the x-independent cutoff A is given by A = 2tcN/L to an accuracy 0(1/ L). The above 
equation determines A A (a;) in terms of A, and clearly AA(x) is x-dependent. For the susy 
kink one can begin by fixing the mode number cut-off for the bosons and A"j ^ for the 
fermions in the trivial sector such that here the bosonic and fermionic mode densities are 
equal. From (JSJ) one then obtains in the kink sector the requirement 



which again determines AA(x) in terms of A. We use this principle to compute the anomalous 
energy density of the bosonic kink, as well as of the supersymmetric kink, which as mentioned 
was obtained already in [Z| through the equality of energy and central charge densities. We 
believe that lmr is sufficient for regularization of the one-loop Casimir energy density in one 
space dimension, and at least necessary in higher dimensions, where further requirements may 
be needed to specify the regularization completely. It is important to emphasize that lmr, like 
ordinary mode regularization, is not easily applied at arbitrary order in perturbation theory, 
because it is not directly expressed in terms of a modification of the action. Thus our claim 
is that for a specialized purpose, namely, calculation of localized Casimir energy, lmr is the 
ideal tool, providing maximal simplicity and efficiency. This claim is simultaneously modest, 
because it is restricted to the computation of one-loop energy densities, and substantial, 
because Casimir energy plays such a large role in quantum physics. 

For the non-anomalous contributions to both the bosonic and susy kink densities, we 
find empirically another striking regularity, phase space factorization. The continuum con- 
tribution to the Casimir energy density in phase space exhibits a remarkable factorization, 
involving a few terms each with simple momentum-dependent factors multiplying functions 
related to the bound-state and zero-mode probability densities in coordinate space. We 
believe this factorization should hold for all reflectionless potentials, but might not extend 
farther. The factorization takes a particularly simple form for what we call below ecas(^), 
a local density whose integral over a region containing the kink gives the total quantum 
correction to the mass of the kink. The local energy density has two contributions besides 
ecas- First, the fundamental definition of Casimir energy density differs from the usual sum 
over zero-point energies by an extra piece which is a perfect differential of an expression 
vanishing far from the kink, so that including this piece does not alter the mass correction 
but does alter the local density. Secondly, there is an effect which perhaps is best viewed 
as vacuum polarization - the quantum fluctuations of the Bose and Fermi fields lead to a 
local shift in the classical background field defining the kink. Again, this changes the local 




(5) 




(6) 
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density but not the total mass. We check explicitly that the energy density and the central 
charge density of the susy kink are equal. 

After these calculations of the anomalous and non-anomalous contributions to the energy 
density of the bosonic (and susy) kink, we turn to the calculation of the central charge density 
of the susy kink. From its definition the central charge density at a point y is an x-integral 
of a bilocal quantity depending on x and y times a delta function S(x — y). Not setting 
x = y too soon yields the anomaly in the central charge density near the kink, confirming 
[7j. Finally, we discuss the physical basis for lmr, observing in particular that point-splitting 
regularization for energy density implies lmr, at least for the bosonic case with arbitrary 
background potential. 

Bosonic kink energy density. — For the energy density of the bosonic kink, one must 
evaluate sums (setting Ti = 1 from now on) ^n^ni^^nix), where the modes <f> n (x) are 
each normalized to unity. As these sums clearly differ from the density sums J2 ( t ) n{ x ) ( t ) n{x), 
one expects in general a nonvanishing one-loop correction to the energy density, and hence 
to the quantum mass. Let us begin with explicit expressions for the mode eigenfunctions, so 
that one may follow the argument in detail. The wave functions of the continuous spectrum 
(using |</> n | 2 (a;) = 1 away from the kink to determine the normalization constant M) obey 



(f)(k,x) 



ikx 



-3 tanh 2 h 1 + 4 

2 



k mx 

+ 6i — tanh 

m 2 



(7) 



with uj = \Jk 2 + m? and J\f 2 = 16^ — ^\ . The zero mode with u = is given by 



<f>o(x) 



1 3m 



V 8 cosh 2 (mx/2) 
The bound state with u B = is given by 



4>b{x) 



'3m sinh(ma;/2) 



' 4 cosh 2 (mx/2) 
The density of the continuous spectrum can be written as follows 



(9) 



\4>(k,x)\ 2 = j^ 



Auj 2 4oj< 



m 2 \ m 2 



-3 



12 



+ 9 



m 2 cosh 2 (mx/2) cosh 4 (mx/2) 
m<f> 2 B (x) 2m0p(x) 



W — 10 



B 



As I 



oo dk 1 
oo 2?r w 2 -u-L 



^ while J™ fH 

m ' J— oo 2-7T u) £ 



— , it is clear that the completeness relation 



2 m 



00 dk 



{|0(fc,x)| 2 -l} + 2 (x)+0|(x) 







(10) 



is satisfied. Eq. (llOj) may be written in a remarkable formula perhaps true for all reflectionless 
potentials, showing factorization of the difference in mode densities in phase space, where the 
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position dependence of each term is given by the corresponding bound-state or zero-mode 
probability density, 

Relation 1: 



2 X m 2 — oj 2 

\<P{k,x)\*-l = -Y.<P 2 M) V ;2 _ / , (12) 

satisfying the completeness relation, as one may check by performing the integration over k. 

Note that all the above expressions for the density do not refer to any particular choice 
of boundary conditions, which of course do affect eigenenergies and the corresponding wave 
functions. The reason is that the choice of boundary conditions will contribute to the density 
away from the boundary at most terms of order 1/L. In the large- L or continuum limit, in 
principle such terms might contribute to the total energy obtained by integration over the 
entire interval between the boundaries. We are unaware of any example of this phenomenon, 
as a previous claim in [2] has been shown to be incorrect 9J. 5 Even if the phenomenon 
were to occur, for the integral just over a finite interval around the kink the effect would be 
negligible, so that the kink energy density and resulting energy can be computed reliably in 
terms of the continuum, modified-plane-wave solutions, unaffected by the choice of boundary 
conditions. 

The requirement in (jSJ) that the topological vacuum density and the trivial vacuum 
density be equal leads via (fTTj) to 



/ — (1 - (k,x )= / dk-y— 2 

J A 71 V ' J A 7T , U) 1 — Ujj J 



71 
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= S(A(l) + ^ (l) ) +0 (^) = ^_^_ + (^) . (13 ) 

With this result for AA(x) we can evaluate the energy density e(x) in the kink sector. 
Adding also the counterterm 6 

I — r A dklm 2 1 3 f A dk , 

AM (x )= Y 2d 2 (x )Jm 2 - u 2 / = ^ — — / — , 14 

y ! Y J° 2nuJ 4 cosh 2 (mx/2) 2vr Jo cu ' V ; 

and rewriting cos ^ mx/2) as ±f B {x) + yields 

6 Cas (x) = e(x) - c W(x) = \ub<&{x) + 2 / A ^(k, x)\ 2 \u - 2 [ A+AA{X) ~u; + AM(x) = 
z Jo 2,71 2, Jo 2tt 2 

5 However, the latter work did identify a delocalized momentum for certain special boundary conditions. 

6 The counterterm AM(x) usually is expressed in terms of the kink background field, but it also can be 
determined by noting that it should cancel the remaining divergence in the difference of the sums over zero 
point energies f_ (dk/2Tv)\<j)(k, ie)| 2 (— j). Equating both expressions yields another formula perhaps valid 
for general reflectionless potentials, 
Relation 2: 

X^ffiaQ^m 2 -^ = 4> 2 K (x)~ ^(oo) 
J2-2jm 2 - uj] JZodxi^ix) -4> 2 K (ao)] 
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1 , . r A dk f md) 2 Jx) 2md£(x)\ AA(i) , / , . 9 , \ f A dk 1 



] -A+mUl(x)+2<f>l(xj) I 

v ' Jo Ztt lu 



2 ray ' Jo 2n\k 2 + m 2 /A k 2 + m 2 J 2vr 

(15) 

The two quadratic divergences proportional to J A dk ui have canceled because we subtracted 
the energy density of the trivial vacuum, while the counter term cancels the remaining 
logarithmic divergence. Again, each term is proportional to a bound-state or zero-mode 
probability density. 

The result is finite and reads 

e Cas {x) = -u B <P B {x)-mJ^ — I - - 1 <f> B {x) - — (4> B (x) + 20 o (x)j . (16) 

The last term is the contribution from the term due to AA(x), and is the analogue of the 
term in the central charge density for the susy case identified as the anomaly by [7]. Using 
the integral 

r°° dk ( 1 1 



we obtain 



o 27r \k 2 + m 2 /4 k 2 + m 2 J 2y/3 



(l m m \ 2 m 2 , 



^ = 7wS > ( 17 ) 



This formula can be rewritten as follows, 




ecas(x) = E i ( 1 - ^ arctan /== % ] "i^O*) - E ^y/rnt-uffiix) , (19) 



where in the first sum the contribution with 1 comes from the bound states, and that with 
arctan comes from the continuum, while the second sum is the anomaly contribution. Such 
formulas for the total mass can be found in [TO] , though we are unaware of local versions 
in the literature. This kink example might be an illustration of a general factorization rule, 
valid for a wide class of reflectionless potentials. While we have not tested it for other cases, 
and do not know how to prove it other than by explicit computation, we believe that its 
simplicity and elegance make the rule worthy of further investigation. 
Integration of ecas(^) over x yields 

M«-> = W-|)-^-^ = ^-^ , (20) 

Z O Z7T 7T 12 ZTT 

in agreement with For convenience later, let us express (|18|) as a total derivative: 



m d 

^Cas 



2 dx I 



< tanh(mx/2) 



9 / / \ ^ 

— tanh (mx/2) 

12 y i > 2 7T 



(21) 



The first term gives the non-anomalous contribution, while the second yields that due to the 
anomaly. 
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Eqs. (J191 l2"T]) give compact expressions for the local energy density, which certainly 
provide the correct total quantum energy of the bosonic kink. However, to obtain the correct 
local Casimir energy density, one must start with an expression for the energy density of 
each mode including a quadratic term in the gradient of the boson field ^(d x rj) 2 , whereas our 
formulae above, giving energies ui n multiplying corresponding mode probability densities, 
implies instead the expression coming from the field equation, — ^r]d 2 r]. 7 Therefore we need 
to add to (fTTJj) the difference, a perfect differential of a function which vanishes far from the 
kink (and so does not change the computed mass correction for the kink), 

Aec-Or) = \d 2 M{x)) . (22) 

The propagator at equal times and positions (t] 2 (x)) (excluding the zero mode which solves 
the homogeneous equation) jTj can be obtained by integrating in (fTUjl w.r.t. * and 

i 2 

adding The divergent part in (rj 2 ) corresponds to setting \<f>(k, x)\ 2 equal to its value at 
\x\ — > oo in the k integral; as the divergence is x-independent it cancels in (|22|) . One thus 

obtains a finite (because the logarithmically divergent part has been subtracted) expectation 
value 

/OO rjh 1 I 
_ oo ^^[l0(^x)| 2 -l] + ^-0|(x) = 

3 I 1 sinh 2 (m3:/2) 



87rcosh (mx/2) 4a/3cos1i (mx/2) 



and therefore, 
d 



Ae Cas = — |tanh(mx/2) 



m/3 1 \ 1 m 



4 V47T 2 y/3 J cosh 4 (mx/2) 16a/3 cosh 2 (mx/2) 



(24) 



As observed in jjj, besides the Casimir energy density there is another consequence of 
the zero-point oscillations, namely, a position-dependent shift in the classical background 
field. This in turn implies a further term in the local energy density, given by 

Ae i4ll) (x) = d x <pid x (j) klnk + (^U 2 )'(f)i = ^(^l^^kink) , (25) 

but of course no shift at this order in the total energy, because the classical energy is sta- 
tionary with respect to arbitrary small variations of the classical field about its equilibrium 
form. Decomposing the Heisenberg field t) as 0kmk(^)+0i(^)+ 7 7(a ; , t), with the quantum 
fluctuation field obeying (rj) = 0, and taking the expectation value of the $ field equation 
(-<9 2 $ + <9 2 $ - (\U 2 )') = gives 

Olh - ( V)"0i = V )"V> - ^Am 2 kmk . (26) 



7 Our boson fluctuation field 77 is designated by x in 0- Note that, as utilized just below, for a mode of 
frequency Lu n our normalized mode function is given by <p n — ^J2u) n rj n . 
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This 0i is just what is needed to satisfy the no-tadpole condition in the kink background. As 
mentioned above, the singularity in (r) 2 (x)) is x- independent and compensated by the Am 2 
term, yielding the quantity (i] 2 (x)) r of (J23j) . 8 Solving (|26|) by the Ansatz 9 0i = v4x0o(x) + 
Bd x <fio(x), and using the fact that 0o is proportional to <9 x 0kmk, one finds 



m 2 




cosh (mx/2) 4 



(md m + 2\d x ] 



-'kink 



(27) 



Through one- loop order (as in the susy case [Zj), the effect of the second term is to replace the 
renormalized mass m and coupling A in 0kink with the bosonic pole mass m = m(l — v^^r) 
given in [Hj, eq.(7) and the adjusted coupling A = ^-A. If we then rewrite the classical 
energy in terms of m and A, the classical mass is multiplied by a factor 1 — y/Sj^-. As 
0i cannot shift the total mass, we know even without explicit calculation that the classical 
energy density in terms of the barred quantities must be renormalized by a compensating 
factor 1 + V3^2. The first term in (|27|) is sensitive only to bosonic fluctuations and hence 
unchanged in the susy case (because, as we shall see, the fermionic source for X includes 
no terms with 1/ cosh 4 (mx/2)); it contributes according to (|23|) . The total one- loop bosonic 
energy density becomes 
Relation 3: 

g ( x ) = [7 2 (A,m,0 kink (A,m,x))^l + v 7 ^— + 

/x a / \ rh / 3 1 \ „ / tanh(mx/2) \ 
+ + A ec „(x) + + —y^ —^-Al ) . (28) 

The last term in (}2"%|) comes from the first term in ([27)1 . The second term in (|27j) renor- 
malizes the classical contribution U 2 , as seen in the first term in ([28)1 . The effect of this 
renormalization is that the classical energy density flattens out a bit. Besides the rescaling, 
all other contributions are of the form d x [<j))^(x)/ cosh™ (mx/2)], with n = 0,2,4. 

Susy kink energy density. — For the susy kink we choose the cut-offs in the trivial sector 
in such a way that the bosonic and fermionic densities in that sector are equal. To make 
the bosonic and fermionic densities also equal in the topological sector, we use a cut-off A 
for the bosons and A + AA(x) for the fermions as in (|HJ). The fermion is described by a 
Majorana two-component spinor ip = (^ + )- As tp + (k,x) is proportional to <p(k,x) while 

ip-(k,x) = ^ (d x + mtanh ip+(k, x) for solutions proportional to exp(— iut) according 
to the Dirac equation jSj, one obtains for the wave functions of the continuous fermionic 
spectrum 

V + (/c,x) = -^0(/c,x) (29) 

8 This procedure gives 4>\ as a finite, renormalized quantity, while \7\ use an unrenormalized <\>\. Hence 
the 'rescaling' part of their <j>\ gives the shift to the pole mass from the unrenormalized mass, but ours gives 
only the shift to the pole mass from the renormalized mass corresponding to the vanishing-tadpole condition 
in the trivial sector. 

9 The term with A is needed for the term in l|23|l proportional to 1/cosh (mx/2) and the term with B is 
needed for terms with 1/ cosh 4 (ma;/2). 
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In the difference of the densities the constant term of course cancels, giving 

\4>(k,x)\ 2 - hMM)l 2 - hMM)l 2 = hMM)| 2 - |V-(M)i 2 

1 ( 9 fuj\ 2 1 \ 

- —727—7^ • (31) 



2J\f 2 \cosh (mx/2) \mj cosh (mx/2) 
For the bosonic and fermionic densities to satisfy © one requires 



+ 2 fl (x) + 2 J o — \4>(k,x)\ 2 



The factor | in |0o(x) comes from the mode expansion i/; + (x,t) = c (j)o(x,t) + ... with 
{co, Co} = 1. 10 This | is the analogue for Majorana fermions of the fractional fermion charge 
discovered by Jackiw and Rebbi for Dirac fermions [TT]. The two terms in parentheses give 
the ip + and ip- contributions of the bound state: tp B+ = \<p 2 B and |?Ab-| 2 = 8cosh ^L x / 2 - ) • We 



obtain 



1 /" A rlk / 

-cj ) 2 {x)+^ B+ (x)-\^\ 2 (x)+2 j( -(|^ + (A,*)| 2 -hMA,*)p 

= (\^ + (k,x)\ 2 + \^(k,x)\ 2 ) . (33) 

Using the completeness relation, and taking the large k limit \ip+(k,x)\ 2 + \i/i-(k, x)\ 2 — > 1, 
one finds 

AA(x) Z" 00 dfc 

7T A 27T 

As we are interested only in the 1/A term, the calculation is easy. From (j31|) we find 



-2 -(\^ + (k,x)\ 2 -\^(k,x)\ 2 ) . (34) 
ie 1/A term, the calcu 
AA(x) m 2 1 



_ (35) 

7r 4nA cosh 2 (mx/2) 

With this result in hand, we compute the difference in energy densities for the susy kink 



+ 2 



ecas,b(x) - e Cas ,f(x) = (fisix) - ip B (x)ip B (x) 
-(|^ + (^x)r-^-(^^)r)^-^^^A + AM susy (x) . (36) 



10 Obtained in |2] from Dirac quantization, this factor equivalently can be deduced from the completeness 
relation for solutions of the single-particle Dirac equation (j>1(x)+2[\^B+\ 2 + \' t pB-\ 2 + J 5§(IV'+| 2 + IV'- 1 2 ~ 1)] = 
0, where \ijj-(k,x)\ 2 = 1/2 — m 2 /(8 cosh 2 (ma;/2)(w 2 — The relative factor of two between the first 

term and the later terms in the completeness relation follows from the fact that, if one sums over a complete 
set of solutions of the Dirac equation, all nonzero frequencies lead to equal contributions from positive and 
from negative frequency, while the zero mode contributes only once. 
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The counter term in the susy case, 



m 2 1 r A dk 1 

2 cosh 2 (mx/ 2) Jo 2n uj 



AM susy (x) = — - I — - , (37) 



is a factor 1/3 smaller, but still nonvanishing. Again, the counter term removes the loga- 
rithmic divergence in the integral, and with (|T7j) one finds 

2 1 m 2 1 

ecas,sus y (V) = e h (x) - e { (x) = (l--)-u B ((p 2 B (x)-ip B (x)ip B (x))-- -j- — , (38) 

3 2 v ' 87T cosh (mx/2) 

where the last term, the contribution from AA, agrees with the central charge density 
anomaly eq.(3.38) in jjj. Integration over x yields the one-loop correction to the mass 
of the susy kink 



M^l y = limx^oo j x dx— 



V3' 



m , , , 3 , ,_ s . , m 



((tanh (mx/2) — tanh(mx/2)) — — tanh(mx/2) 



48 u v ' J y ' " 4tt 

1 . . m m , , 

= -u B {l-l - _ = __ , (39) 

which of course is the accepted answer. Note that the non-anomalous contributions from the 
bosons and the fermions do not cancel locally, but in the integral they do: |^b(1 — 1) = 0. 
For explicit expressions later it is helpful to rewrite the first part in the bracket of (J39|) as 

d TTt 

ecas,sus y (non — anom) = 1=(~ tanh(mx/2)/cosh 2 (ma;/2)) . (40) 

dx 16v3 

As in the bosonic case we must add the missing term in the bosonic Casimir energy 
density Aec as given in eqs. ()22I23|) . as well as include the shift for the susy case <pi(x) in 
the background field. We compute this 0i, again using the second-order field equation for 
$, which now has an additional contribution from fermions: 



dl<Pi ~ (^ 2 )>i = {^(\ U T(V 2 } ~ ^Am 2 kink } + {lu"(W) - iAmj0 kink } , (41) 
with 



[-1 + 6m 2 /cosh 2 (ma;/2)]tanh(ma;/2)/16(A; 2 + m 2 /4)] 

-do 2ttuj 

- v / 3m/4tanh(mx/2)/ cosh 2 (mx/2) . (42) 

The last term comes from the bound state, and the term with (k 2 + m 2 /4) in the numerator 
of the integrals is cancelled by the fermionic part of the mass counter term. Again using 
(117)1 . we find that the fermionic contributions to 4>i are only proportional to 1/ cosh (mx/2), 
and not 1/ cosh (mx/2), so that the net coefficient of the 1/ cosh 2 (mx/2) term is a factor 
2/3 smaller than in the bosonic case. The final result for <pi reads 



m z 



--(md m + 2\dx 



2^3 47r/ cosh 2 (mx/2) 2^3 



(43) 
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Note that in (}4*3*j) the first term is the same as in the bosonic case (as mentioned earlier), 
while the second term is smaller by a factor 2/3. This is the same result found in |7] using 
a first order differential equation based on susy considerations. Iterating the susy relation 
(d x (p + U) — 0, one confirms that the second-order and first-order approaches are consistent 
with each other. 11 For the susy energy density we then find full agreement with the central 
charge density of [7j, after restoring a missing factor of \ in the first line of (5.21) in that 
work, kindly pointed out to us by the authors. 
Relation 4: 



£(x) = £/ 2 (A s ,m s ,0 kink (A s ,m s ,x))(l + 



A,, 



ml 



, s m s / 3 1 \ _ / tanh(m s x/2) \ . . 

+ e<w*) + a^m + T ( s + ^H C0Sh \ m J^ ) ■ m 

Here we use m s = m(l — \/2\/3m 2 ) and A s = A(l — \/\^3m 2 ). Thus the rescaling of the 
classical density involves a shift 2/3 as big as for the bosonic case, while ec a s,sus y includes 
an anomaly 1/3 as big as in the bosonic case and a different finite-energy contribution, as 
explained above. Meanwhile, the last two terms, from Aecas and from the nonrescaling term 
in 0i r , are the same as in the bosonic case. 12 

Central charge density. — We now compute the anomalous contribution to the density 
((x). Before regularization one has H(x) = \4> 2 + \4>' 2 + \U 2 + ^ (V'+V'h- + V'-VO ~ iU'i/) + xp^, 
and ((x) = U d x (j> (note that [7; have the opposite sign convention for (). Using the equal-time 
anticommutators of the fermionic fields ip+{x) and ip-(x) and the definition 



J± 

one obtains 



+ C/hV± , (45) 



{Q±, j±(y)} = 2H(y) ± 2({y); j ± = + (<P' ± U)^ T ; Q ± = ]± dx ; (46) 

J-L/2 



C(v) = / dx 



\ {{^ + {x)^ + {y)} - {^_(rc), V-(y)}) (\kx)ky) ~ y^Wiv) + \u{x)U{y) 



(47) 



+ l -{{i, + {x)^ + {y)} + {i,.{x)^{y)}) fy{x)U{y) + hj{x)<t>\y) 



n The details are as follows. From (d x c/> + U) = we have d x <\> x + U'fa + (1/2!)C/"(t ? 2 ) + v /A/2(-A A i 2 /A) = 
. Differentiating with respect to x, using c^kink = — U, and eliminating d x 4>\ yields an equation for 9 2 0i. 
That this is equivalent to l|41() follows from the identity (rj(d x + U')r] + ((d x + U')rj)rj) = — 
which in turn is a consequence of ip- — (i/uj)(d x t&nh.('mx/2))ip + and -0+ = vV^- 

12 For comparison with note that our coupling A is equal to 2A 2 in their formulation. Also, they give 
the integral of the density from —x to x, while we write the pieces of the quantum correction to the density 
as local derivatives, so that our expression for the function being differentiated is half theirs for the integral. 
The specific terms in their equation (5.21) are related to ours as follows: The first line in (5.21) is simply the 
integral of what we call U 2 . The first term in the second line is the anomaly. In the final bracket, the first 
term receives equal contributions from ec as and Aec as - The remaining piece receives equal contributions 
from Aec as and from the nonrescaling part of (f>i, the shift in the classical field. 
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This is an exact result; all terms with bosonic commutators cancel. The anticommutators 
{if) + (x),if)-(y)} and {ip + (y),ip-( x )} all vanish, and also the first line in ((y) vanishes, while 
the second line would seem to give 



/+oo 
dx5(x -y) 
-oo 



-(r](x)r](x))U"<p'(x) + (r)' \x)r]{x))U' \x) 



d 

dx5(x-y) — 



-(r)(x)r)(x))U'(x) 



(48) 



DO ] 

dy((y) = 7;(v(x)v(x))U'(x)\™ 00 . (49) 

oo / 

This is the expression obtained in [3]. Below we show that with appropriate care (i.e., not 
setting x = y too soon), there is an extra term - the anomaly. The naive result in (}4*9*|) 
contains a free field propagator for rj, because at x = ±oo the effects of the kink disappear, 
and, adding the counterterm to the central charge due to mass renormalization, all quantum 
corrections to the central charge would seem to vanish. In the approach of 7J, on the other 
hand, the central charge contains a naive term (p'U and an explicit correction term which 
is also a total derivative and proportional to 1/M 2 . Because their rj propagator contains 
an extra regulating factor (k 2 + M 2 ) -1 , the contribution in (|49J) now cancels even after 
regularization, but because the correction term contains two extra derivatives (to balance the 
factor 1/M 2 ) it yields an extra contribution proportional to M 2 /M 2 , which is the anomaly. 

In our case we start from (|4T|) . but without extra terms as in We keep x ^ y in (}4Tj) . 
giving 

r 11 
C(y) = J d x5(x-y)[U'\x)^ kink (y)-(r] 2 (x)) + -(ri'(y)ri(x))U'(x) + (x ?/)+A/i 2 term] . (50) 

We now show that the result is still a total derivative, but instead of the total derivative 
in (@HJ), rather a total derivative with an extra term. The crucial point is that one cannot 
replace 5(x — y) (i]' (x)rj(y)) by \5{x — y)d x {rj 2 (x)) because there is a singularity in (t)' (x)r](y)) 
as x tends to y. Setting 5(x — y)[{i]'(x)ri(y)) — (l/2)d x (r](x) 2 ) = would mean that all terms 
vanish as in j^J. However, the singularity which invalidates this equality yields the anomaly 



((x)dx — 



47T 



(51) 



where W'{<j>) = U. Hence M« = -Z^ in agreement with the invariance of the background 
under Q + (which corresponds to the susy transformation with parameter e_). Let us see 
this explicitly. 

The identity we need is 

S{x - y)( V '(x) V (y))(f(x) - f(y)) = -±-6{z - y)f(x) , (52) 

where / is any smooth function of x. The proof of this identity follows from (7](x)r](y)) = 
— 7^ In \x — y\ + A(x, y), where A is a smooth, symmetric function, therefore involving only 
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nonnegative even powers of (x — y), as can be seen from (J7EJ). The actual calculation of 
the anomaly is now very simple. Expanding all contributions in terms of x — y, and using 
S(x — y)(x — y)Afi 2 = after regularization of A/z 2 , while 5(x — y)(x — y)d x (rf(x)) = does 
not need regularization because d x {rf(x)) is finite, we have from (JH2|) . 



/OO rOO 
C(y)dy= / dxdy 
-oo J—oo 



^5{x-y)U"{4>W{^) 

47T 



s ■ (53) 



Here we used C/ = y| (</> 2 — ^-j, = ^=tanh^|, and m = \i\f2. Again we have the 
accepted result, and we see that point-splitting regularization yields the same extra term in 
the central charge as does higher-derivative regularization. 

Thus we have shown in a simple way that the term (rj'(y)rj(x))U'(x) produces the anomaly 
if one does not set x = y too soon. In [7j a more complicated but also more powerful 
regularization scheme was used to prove this. Our observation pinpoints the place where 
naive methods missed the anomaly. As discussed extensively in the previous sections, one 
must add to the anomalous part the various contributions to the non-anomalous part of the 
central charge density. This works exactly as in 7j, and of course is completely unaffected 
by choice of regularization method. 

In view of our emphasis on lmr for energy density, it is reasonable to ask why we do not 
attempt to apply it to central charge density. Looking at (jEDJ), one sees that the expression 
to be regulated, the bilocal correlator in i](x) and i]'(y), which clearly is not determined 
by insisting that the regulated sum lm{rj(x)rj(x)) is unchanged between vacuum and kink 
backgrounds (the entire content of the lmr prescription). Thus lmr may be applied as a 
condition on the expressions in the central charge density, but is not sufficient to regulate 
them. 

Foundations and conclusions. — Finally we comment on the physical basis for lmr. In 
Planck's original formulation of quantum physics ^2], the number of degrees of freedom is 
defined by the available volume in phase space. To fix the total number of modes while intro- 
ducing a background potential affecting the fluctuations is simply to conserve the total phase 
space available. The work of Einstein ^3] and Debye [Hj on crystal vibration contributions 
to heat capacity introduced the concept of a local density of degrees of freedom, codifying a 
notion already found in Boltzmann's lectures on gas theory . As was true for their work, 
in a lattice approach the number of degrees of freedom per unit volume evidently does not 
change when interactions are introduced, and the local mode density should be equal to this 
number of degrees of freedom. 

Also point splitting methods clarify the meaning of lmr. Consider the bosonic local mode 
density regulated by point splitting 

p(k, x) = |" dy<f>* (k, x - |) <j> (k, x + |) f(y) , (54) 

where f(y) is a function sharply peaked around y = 0, with J dyf(y) = 1. For large k, the 
JWKB approximation for <ft(k, x) is 

4>{k,x) =e ^ e -iJ x dx'v( x ')/2k ^ ( 55 ) 
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where V(x) = U(<f>(x))U"(<p(x)) + (U') 2 ((f)(x)) - (U'f((/)(\x\ -> oo)). Substituting this ex- 
pression into p(k, x) one finds for the integrand of 



e iky e 2k f(y) ~ e< k —^) v f{y) . (56) 

In the trivial sector p(k,x) = f(k), where f(k) is the Fourier transform of f(x), but in the 
kink sector one finds a modification p(k, x) = f (k — V(x) /2k). The energy density therefore 
contains a term 8e(x) = J 5p(k, x)^u>7jr, and expanding / we find 13 

For the bosonic kink, this is the anomaly in ()16|) . as one may readily check by direct substi- 
tution. 

From the JWKB form for the wave function at high energies it follows that the quantity 
AA(i), and hence the anomaly, depends on x only through the potential felt by the fluctu- 
ations. This in turn implies that the local anomaly in the energy density is determined at 
each x by the value of the classical background field <p(x), as stated for the central charge 
density of the supersymmetric case in [Jj. 

While the above discussion shows that the lmr result for the bosonic kink follows from 
point-splitting, it is possible to make a much stronger statement, that point-splitting implies 
lmr for a scalar field in an arbitrary background potential in one space dimension. The 
same JWKB approximation used to determine the shift in effective wave number due to 
the potential V can be used also to compute the modulation of the mode density versus 
energy (or equivalently, versus asymptotic wavenumber k) in the region of nonzero potential, 
compared to the asymptotic density far away. This is a standard calculation, obtaining the 
wave function and hence the density to one higher order in 1/k than required for the shift in 
effective wavenumber. There is a simple physical mnemonic for the result of the calculation. 
Treating this system as a Schrodinger problem with "Hamiltonian" u 2 , the modulating factor 
is simply the ratio of the asymptotic classical velocity to the local velocity: 

p(x)/p(oo) = v/v{x) = k/^Jk 2 - V(x) , (58) 



or 

5p/p ~ V(x)/2k 2 . (59) 

Integrating the above expression from a nominal sharp cutoff A to oo, we see that to have 
the same integrated mode number density above the cutoff in the presence of a potential 
V(x) as in a trivial background, we must shift the cutoff by 

5A = V(x)/2A , (60) 

13 The narrower f(x), the wider the range in k contributing to (|57J) . so that in the limit the contribution of 
any finite range around k = becomes negligible. This justifies our use of the JWKB approximation, which 
is valid for k 2 ^> m 2 , V. 
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exactly the amount implied by point-splitting as found in the discussion leading to (|57jl 
above. 

The equality between the wave number shift 5k(x) and the integral over the density 
shift is reminiscent of an unsubtracted dispersion relation. Possibly the implementation of 
lmr in higher dimensions would require the equivalent of subtracted dispersion relations to 
compensate for the increasing divergence of energy density with cutoff. 

We have seen that point-splitting, a regularization scheme which is local but not neces- 
sarily useful beyond one-loop order, gives the same anomaly in the central charge found in 
[Z| with higher-derivative regularization, and also implies lmr for the bosonic energy density 
with arbitrary background potential. As explained above, the converse is not true: lmr does 
not contain the full content of point-splitting regularization. Nevertheless, it is appealing 
that it captures in one line the above sequence of equations required to obtain the anomaly 
from point-splitting regularization. Thus, as stated already in the introduction, lmr is a 
simple, easily-used tool for a special but extremely important application, the computation 
of local one-loop energy densities. 

For the supersymmetric case, one gains insight into the requirement of equal bosonic and 
fermionic mode densities by considering the N = 2 theory, where there is an abelian charge 
density which should be invariant under supersymmetry. That requirement automatically 
imposes the constraint represented by lmr. 

While all of the above are appealing arguments, the accepted criterion for determining 
the validity of a regulation procedure is to insert regulators into the action in such a way that 
all relevant symmetries are satisfied at the regulator level, and then to deduce consequences 
for specific quantities. Thus in the present case a definitive check on the validity of lmr 
would be to use, for example, higher- derivative regulation (which obeys supersymmetry), 
and check that this scheme implies lmr. This important analysis remains to be done. As an 
alternative, one might be able to prove that point-splitting preserves susy in models with 
solitons, and then extend our deduction of lmr from point-splitting in the bosonic case to 
the susy case. This is something to which we intend to return in the future. 

To summarize, lmr permits one to isolate and then compute directly the anomalous con- 
tribution to the energy density of the bosonic or susy kink. Expressed most conservatively, 
lmr at the least gives a simple interpretation of the anomaly as the shift in energy density 
required to equalize mode densities. In fact as we have just seen, at least for Bose fields in one 
space dimension with arbitrary scalar background potential, lmr follows from point-splitting 
regularization of the energy density. In addition we found remarkable phase space factor- 
ization identities for the non-anomalous contributions to the energy density, which might 
hold for all reflectionless potentials. These non-anomalous contributions are independent of 
the regularization method (though sensitive to renormalization conditions because they only 
are convergent after subtraction of the mass counter term). Elsewhere [H>| we compute the 
divergent energy density at the boundary of the kink with supersymmetric boundary con- 
ditions, and obtain an analytic expression for the anomaly near the boundary, which in the 
limit when the regulator energy goes to infinity becomes a delta-function contribution just at 
the boundary, in agreement with expectations of [7j. It would be interesting to explore local 
mode regularization further, comparing with complete regularization schemes and studying 
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solitons in higher dimensions such as the magnetic monopole, and also explore phase space 
factorization, seeking a theoretical basis as well as additional examples. 

We thank Kevin Cahill, Ludwig Faddeev, and Stefan Vandoren for discussions, and in 
particular Arkady Vainshtein for useful detailed comments on a number of occasions. We note 
an independent proposal of lmr by Robert Wimmer which appeared after submission 
of the present work. 
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